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I.  Introduction 


This  is  the  AASERT  Final  Technical  Report  on  work  performed  under  the 
support  of  the  Air  Force  Office  of  Scientific  Research  under  Grant  F49620-93-1- 
0366  for  the  period  June  1,  1993  to  November  30,  1996.  The  objective  of  this 
work  was  to  train  graduate  student  and  teach  him  how  to  use  particle-in-cell 
code  to  simulate  the  process  of  radiation  generation  from  relativistic  electron 
beams.  At  the  same  time,  we  have  developed  a  new  version  of  particle-in-cell 
code. 

2.  Summary  of  Work  Accomplished 

In  this  work,  we  have  developed  a  new  version  of  particle-in-cell  code 
which  is  capable  of  separately  following  the  evolution  of  electromagnetic  and 
electrostatic  fields.  The  new  code  was  utilized  to  investigate  the  space  charge 
effects  in  gyro  backward  wave  oscillators. 

The  electron-wave  interaction  in  microwave  devices  is  a  highly  nonlinear 
phenomenon  and  reguires  computer  simulations  to  follow  its  evolution.  The 
electron  convection  current  calculated  in  simulation  often  contains  both  the 
electromagnetic  and  electrostatic  components.  This  total  current  is  usually  used 
as  the  source  for  exciting  the  electric  and  magnetic  fields  in  electrodynamic 
structures.  In  this  approach,  it  is  important  to  make  sure  that  the  equation  of 
continuity  is  satisfied  at  all  time.  In  MAGIC\  this  is  accomplished  by  correcting 
the  electric  field  at  every  time  step  to  insure  that  Poisson's  equation  is  satisfied. 
Most  of  the  other  approaches  employed  to  simulate  microwave  devices  use 
medium  electron  beam  current  and  assume  fixed  transverse  field  profiles 
determined  from  the  boundary  conditions  of  the  structure.  The  electrostatic  field 
evaluated  from  the  Poisson  equation  is  often  added  to  the  fields  generated  by 
the  total  current.  In  doing  so,  the  dominant  mode  of  space  charge  fields  is 
doubly  counted^.  In  order  to  be  able  to  correctly  investigate  the  effects  of  space 
charge  field  on  the  performance  of  microwaves  devices,  we  have  developed  a 
new  version  of  particle-in-cell  code  which  is  capable  of  isolating  the  transverse 
current  from  the  total  current.  The  new  code  is  used  to  simulate  gyro  backward 
wave  oscillators  and  the  result  is  presented  in  the  report.  In  the  future,  the  same 
algorithm  will  be  applied  to  simulate  plasma-filled  backward  wave  oscillators 
and  helix  traveling  wave  amplifiers. 
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A.  Algorithm  for  Isolating  Transverse  Current 

By  using  the  Coulomb  gauge^,  the  current  density  (or  electric  field)  can 
be  written  as  the  sum  of  two  terms, 

+  .  (1) 


Ji  is  called  the  longitudinal  (irrotational,  V  xJf  =0)  current  and  is  the  driving 

source  for  electrostatic  fields,  while  jt  is  called  the  transverse  (solenoidal, 

V«jf=  0  )  current  and  is  the  driving  source  for  electromagnetic  fields.  The 
magnetic  field,  of  course,  only  has  transverse  components.  Starting  from  the 
vector  identity, 

Vx(Vx‘?)  =  V(V.'?)-V2^  ,  (2) 

the  transverse  current  can  be  determined  from  the  total  current 

v2(‘?,)=-VxVx't  =  v2Tf-V.(  V.Tf  )  (3) 


The  solution  of  Eq.  (3)  in  integral  form  is^ 


^t=— VxVx 
4k 


't') 


d3x' 


(4) 


This  form  of  solution  is  difficult  to  solve  numerically.  In  the  differential 
formulation,  Eq.  (3)  can  be  easily  solved  if  the  field  quantity  is  expanded  into  a 
series  of  complete  orthogonal  functions.  For  example,  sets  of  such  functions 
are  the  Fourier  series  and  the  Bessel  functions  which  are  the  solutions  for 
cylindrical  waveguides.  Consider  an  infinite  and  periodic  system.  In  this  case  it 
is  convenient  to  solve  Maxwell's  equations  in  the  Fourier  space‘s  and  the 
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spliting  of  the  current  into  transverse  and  longitudinal  components  can  be  easily 
done.  From  Eq.  (3),  the  transverse  current  can  be  written  as 


-f  -f 


(5) 


where  k  is  the  wavevector.  Based  on  this  approach,  for  the  last  two  decades, 
we  have  unravelled  numerous  interesting  nonlinear  phenomena  in  plasma 
physics. 

In  microwave  devices,  a  cylindrical  waveguide  (r,6,z)  is  often  provided  to 
confine  the  electromagnetic  wave  propagation  along  the  axial  direction.  A 
combination  of  Bessel's  functions  can  always  be  found  to  satisfy  the  proper 
boundary  conditions  and  Maxwell's  equations.  We  shall  first  develop  an 
algorithm  under  the  assumption  of  axial  symmetry.  In  this  situation,  the 
electromagnetic  field  profiles  can  be  expressed  in  terms  of  the  TEon  and  TMon 
modes.  Since  the  electric  field  in  the  TEon  mode®  is  only  along  the  0  direction 

and  is  independent  of  the  0  coordinate  (V»~^=0),  it  contains  only  the 
electromagnetic  component.  For  the  TMon  mode,  the  electric  fields  are  along 
the  z  and  r  directions.  To  obtain  the  transverse  current  from  the  total  current,  Eq. 
(3)  can  be  used.  Take  the  spatial  variation  of  the  electron  beam  current  to 
conform  the  electric  field  variation 

^  =  [zJzJo(Pon^)+  r  JrJl  (  Pon  ^  )]  (6) 

where  a  is  the  waveguide  radius  and  pon  is  the  zero  of  the  Bessel  functions  Jq. 
The  first  term  of  the  right-hand  side  of  Eq.  (3)  becomes 


The  transverse  current  can  be  written  as 
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Ju  =  -  ( #  Jz  +  ^  J, )  Jo(  Po„  ^ )  . 


a" 

2 


(8) 


J,,  =  -(kz^J,-^Jz)Ji(Po„i)e*^^/(^  +  kz=’) 


Notice  that  Eq.(8)  is  divergence  free.  These  currents  are  retained  as  the  source 
to  advance  the  electromagnetic  fields  through  Maxwell's  equation. 


at 


=  c  V  X  #  -  47r  jt 


(9) 


-  =  -cVxE. 


In  the  next  section,  the  electrostatic  field  will  be  evaluated  from  Poisson's 
equation. 


B.  Using  B-Spline  to  Solve  Poisson's  Equation 

For  electromagnetic  wave  problems,  the  field  in  the  exterior  of  the 
electron  beam  takes  the  form  of  a  wave  moving  outwards  if  there  is  no  structure 
to  confine  the  wave  propagation.  In  a  cyclindrial  waveguide,  the  transverse 
mode  is  a  global  mode  and  one  mode  is  often  sufficient  to  represent  the 
electron-electromagnetic  wave  interaction.  For  electrostatic  problems  the  fields 
outside  the  beam  location  decay  exponentially,  with  no  oscillation  spatially  and 
no  wave  propagation.  To  simulate  this  localized  field,  a  great  many  transverse 
modes  are  needed  along  the  radial  direction.  The  scheme  of  using  finite- 
differencing  or  B-spline  may  be  a  better  approach  to  evaluate  the  electrostatic 
field.  Here,  we  shall  describe  the  B-spline  scheme  for  solving  the  Poisson 
equation. 

The  two-dimensional  Poisson  equation  in  cylindrical  coordinates  (r,z)  is 


8^<l>  1  3<I>  3^0 

3r2  ^  3r 


=  -  4jT:p  . 


(10) 
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The  boundary  conditions  are  O(r,0)  =  ^>(r.L)  =  <E>(a,z)  =  0.  Expand  the  functions 
O  and  p  along  the  z-direction  in  terms  of  sine  function? 


(11) 


where  Nz  is  the  number  of  grid  points  in  the  z-direction.  Substituting  Eq.(11) 
into  Eq.(10)  yields 


^k"  +  7  <E>k'  -  (^  )^  Ok  =  -  47rpk  ,  (12) 

where  prime  denotes  taking  derivative  with  respect  to  the  r  coordinate. 

The  collocation  method  with  B-splines®  is  used  to  numerically  solve 
Eq.(12).  Let  {  n  be  the  radial  grid  points  (knots),  0  =  ri  <  ....<  r^^  =  a.  The 

solution  of  Eq.(12)  can  be  approximated  by  a  spline  S  of  fourth  order  on  the  grid 
{  n  }.  A  spline  can  be  expressed  in  the  following  form 


S(  r)  =  ^aiBi(r) 


i=1 


(13) 


where  Bi(  r )  are  so  called  B-splines  which  are  basic  functions  of  the  space  of  all 
cubic  splines  on  the  grid  {  n  }.  n  =  Nr  +  2  is  the  dimension  of  this  spaced. 
Substituting  S(  r )  for  <I>(  r )  in  Eq.(12)  we  obtain 


i=1 


(14) 


Since  the  values  of  p,^  are  known  at  the  grid  point  q,  we  will  evaluate  Eq.(14)  at 

(k) 

these  points  to  determine  the  coefficients  af  . 
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^  / 1-\  II  I  2 

X  ar^[Bi  (rj)+-Bi(rj)-(-)  Bi(  rj )  ]  = -47tpk(  n )  , 

i=i  ^ 

j  =  2,3, . Nr  (15) 

<I>' 

Because  the  term  —  in  Eq.(IO)  is  singular  at  r  =  0,  we  must  treat  it  carefully.  The 

property  of  axis-symmetry  implies  Er  =  -O'  =  0  at  r  =  0.  The  L'Hospital  rule  can 
be  used  to  evalute  this  singular  term 

Iim7  0'=  lim  O"  =  O"(0)  .  (16) 

r-40  '  r^O 


Using  Eq.(16),  Eq.(14)  at  r  =  0  becomes 

n  /ix'X  1 1  2 

£  ar'[2B|  (O)-(-j-)  Bi(0)l  =  -()iPk(0)  .  (17) 

Besides  Eq.(15)  and  (17),  we  have  two  additional  boundary  conditions  <I>'(0)  =  0 

n  ...  n 

and  Ok(  a )  =  0.  In  terms  of  B-splines,  these  become  ^  ai' '  Bi'(  0  )  =  0  and  ^ 

i=1  i=1 

Bi(  a  )  =  0.  Now  we  have  altogether  Nr  +  2  equations  for  Nr  +  2  unknown 
coefficients  af  and  the  resulting  equation  can  be  written  in  the  form  of 
tridiagonal  matrix.  After  have  been  evaluated,  the  electrostatic  field  can  be 

determined  from  =  -V  x  O. 


The  new  version  of  particle-in-cell  code  described  in  the  previous 
sections  was  used  to  simulate  the  performance  of  gyro  backward  wave 
oscillators.  The  parameters  (table  I)  used  in  simulations  is  also  roughly  the 
parameters  employed  in  a  recent  experiment  carried  out  at  Phillip's  Laboratory. 
The  experimental  results  show  very  low  efficiency  (less  than  0.1%)  when  the 
TEoi  mode  was  used. 
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Table  I:  Parameters  of  the  gyro  backward  wave  oscillator 


Waveguide  radius  (rw) 

4.37  cm 

Operating  mode 

TE01 

Operating  cyclotron  mode(s) 

1 

Beam  voltage  (Vb) 

400  kV 

Beam  current  (Ib) 

2kA 

Guiding  center  position 

0.48  rw 

Electron  velocity  ratio  a  (vi/vz) 

0.7 

Magnetic  field  (Bo) 

3.7  kG 

System  length  (L) 

20  cm 

Figure  1a  displays  the  dispersion  relations  of  the  TEoi  mode  and  beam 
cyclotron  lines.  The  intersection  point  "A"  is  the  desired  backward  wave 
oscillation.  Using  a  cold  and  zero  spatial  spread  beam,  the  beam  annular  width 
is  only  twice  the  beam  Larmor  radius.  The  beam  plasma  frequency  (cob)  can  be 
evaluated  and  is  about  0.524  of  coc  (waveguide  cutoff  frequency).  Figure  1b 
shows  the  time  evolution  of  output  power  determined  from  simulations.  Without 
taking  into  account  the  effects  of  space  charge  field,  the  saturation  efficiency  is 
about  2%.  Under  the  influence  of  the  space  charge  field  arising  from  the 
neutralized  beam,  the  output  power  reduces  slightly  while  the  growth  rate 
remains  relatively  unchanged.  This  is  because  in  a  gyrotron,  the  azimuthal 
phase  bunching  dominates  over  spatial  bunchings.  The  AC  space  charge  field 
plays  a  minor  role.  Figure  2  plots  the  spatial  distribution  of  the  beam 
electrostatic  potential  and  electric  fields.  Outside  the  beam  region,  the  field 
quantity  does  decay  exponentially.  The  experimental  result  from  Phillip's 
Laboratory  indicates  that  the  efficiency  is  much  better  (almost  10  times)  if  the 
TE21  mode  is  the  operating  mode.  In  order  to  simulate  this  results,  a  three 
dimensional  electrostatic  code  is  required.  We  intend  to  develop  such  a  code 
by  including  a  finite  number  of  azimuthal  modes.  We  are  also  working  on 
electrostatic  codes  to  simulate  plasma-filled  backward  wave  oscillators  and 
helix  traveling  wave  amplifiers.  Both  devices  require  periodic  waveguide  wall 
to  slow  down  the  electromagnetic  wave  propagation.  We  will  also  investigate 
the  effects  of  a  non-neutralized  beam  on  the  performance  of  gyro  backward 
wave  oscillators. 
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Figure  2  Transverse  spatial  distribution  of  electrostatic  field  (k=1),  (a)  potential, 
(b)  radial  electric  field  (Er),  and  (c)  longitudinal  electric  field  (Ez). 


